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ABSTRACT 



A large number of carbon-rich evolved objects (asymptotic giant branch stars, 
protoplanetary nebulae, and planetary nebulae) in both the Milky Way galaxy 
and the Magellanic Clouds exhibit an enigmatic broad emission feature at ~ 
30 /mi. This feature, extending from ~ 24 fim to ~ 45 /im, is very strong and 
accounts for up to ~30% of the total infrared luminosity of the object. In 
literature it is tentatively attributed to magnesium sulfide (MgS) dust. Using 
the prototypical protoplanetary nebula around HD 56126 for illustrative purpose, 
however, in this work we show that in order for MgS to be responsible for the 
30 jttm feature, one would require an amount of MgS mass substantially exceeding 
what would be available in this source. We therefore argue that MgS is unlikely 
the carrier of the 30 /zm feature seen in this source and in other sources as well. 



Subject headings: dust, extinction - circumstellar matter - infrared: stars - stars: 
AGB and post-AGB - stars: individual (HD 56126) 



Introduction 



In carbon-rich evolved objects, there are two pr ominent, mysterio us emission features 
known as the "21 fim" and "30 /zm" features (e.g. see I Jiang et al.ll2009l ). The 21 /im feature 
is seen almost exclusivel y in 16 protoplan e tary nebulae (PPN e) and its identification is 
notoriously difficult (see iPosch et al.l 12004 : IZhang et al.l 120091 ). The 30 /im feature, first 
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discovered by iForrest et al.l (Il98ll ) in the Kuiper Airborne Observatory spectrometry of C 
stars and planetary nebulae (PNe), is very broad and strong, extending from ~24/im to 
~ 45 Atm and acco unting for up to ~ 30% of the total infrared (IR) luminosity of the object 
flVolk et alJliooi lfl More ubiquitously seen in C-rich objects than the 21 fim feature, the 



30 /im feature has been detected in 63 Galactic and 25 Magellanic objects (see iJiang et al. 



20091 and references therein), including asymptotic giant branch (AGB) stars, post-AGB 
stars and PNeEl 

Compared to the 21 /im feature for which over a dozen carrier candidates have been 
proposed, the identification of the 30 /mi feature has not received the attention it deserves. 
Magnesium sulfide (MgS) solids were first proposed by lGoebel fc Moseleyl ( 119851 ) as a carrier 
of the 30 /im feature, based on (1) the similarity of t he observed emission spectral profiles 
and the laboratory spectra of MgS (INuth et al.lll985l ). and (2) the considerations of C-rich 
equilibrium condensa tion chemistry which predicts the condensation of MgS in C-rich stars 
(ILattimer et al.lll978l ). The MgS proposition gains further support from detailed modeling 
of the 30 /tm featu re of a large number of C-rich sources obtained by the Infrared Space 
Observatory (ISO) j,Tiang et alill999l : ISzczerba et alJEoOflL Irlonv et alJl2002h . 



A valid carrier candidate should not only reproduce the observed emission spectra but 
also satisfy the abundance constraints (i.e., the amounts of elements required to lock up 
in the carrier should not exceed what are available in the 30 /im sources). The absorption 
profile Qabs(A) of MgS dust with a distribution of ellipsoidal shapes exhibits a broad band 
around 30 /im. To fit the observed 30 /im emission spectra, one often compares the observed 
spectra with Qabs(A) B\(T), the pro duct of the absor ption efficiency and a blackbody with 
an assumed temperature T (e.g. see lHony et al.l 120021 ) . 



The MgS temperature is an essential parameter: it not only affects the shape and 
peak wavelength of the model emission spectrum but also determines the total amount 
of MgS dust MMgs required to account for the flux Fx emitted from the 30 /im band: 
Myigs oc F\/ [Qabs(A ) B\(T)}. Previous models based on MgS do not seem t o have a sulfur 
(S) budget problem (I Jiang et al.lll999l ; ISzczerba et al.lll999l ; iHony et al.ll2003l ). However, the 
MgS temperatures adopted in these models were often treated as a free parameter and were 



1 The ISO spectroscopy suggested that the 30 /im feature consists of tw o subfeatures: a narrow feature 
at 26 /i m and a broad one a t 33 /im (jrlrivnak et alJl200(]Hvolk et alj|2002h . But this was not confirmed by 
Spitzer (jflrivnak et al.ll2009h . 

2 The 21 /im sour ces all exhibit the 30 /im feature. But these two features do not appear to correlate in 
their strengths (see Jiang et al. 20091 ). Unlike the 21 /im feature which displays little shape variation, the 
30 /im feature var ies in its peak wavelength and width among different sources (e.g. see Hrivnak et al. 2000t 
Honv et alj|2002h . 
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derived from matching the spectral profil es of the observed 30 fim feature with Q a b s (A) B\(T) 



by varying T (e.g. see iHony et al.l 120021 ) . Due to the lack of the dielectric functions in the 



ultraviolet (UV), visible and near-IR wavelength ranges for MgS dust, it is not possible to 
calculate the precise thermal equilibrium temperatures of MgS dust in the 30 /im sources. 
Therefore, the previously derived S budget may be questionable. 

In this work we make a rather generous estimation of the UV/visible absorptivity of 
MgS and calculate its temperatures in the circumstellar shell of HD 56126, a prototypical 
source of the 21 /im and 30 /im features. It is shown that the required MgS dust mass exceeds 
what would be available by a substantial amount. We therefore argue that MgS is unlikely 
the carrier of the 30 /im feature. 



2. Optical Properties of MgS Dust 

In order to calculate the temperature of a dust species, the knowledge of its optical 
properties over a wide wavelength range is required. Unfortunately, the di electric functions 



of MgS have only been exp erimentally determined in the IR (e.g. see iBegemann et al. 



ixp* 

19941 ; iHofmeister et all 12003 ). The UV/visible absorptivities of MgS dust are essential since 



they determine how much energy will be absorbed by a given MgS grain (when exposed to 
starlight) which in turn determines how much energy will be emitted in the IR. It is this 
balance of energy between absorption and emission that determines the grain temperature. 



We take the following procedure to approach the UV/visible absorption efficiency Qabs(A) 
of MgS dust: 

gS(A) = { to (Ao/A) ^ "A > A ; (1) 
where Q a is a constant and the cut-off wavelength A Q depends on grain size (and material). 
What would be the reasonable Q D and A Q values for MgS? To shed light on this, we use 
Mie theory to calculate the UV to near-IR absorption efficiencies for spherical grains of 
various compositions and sizes. We consider both dielectric (amorphous silicate, amorphous 
carbon, organic refractory) and metallic materials (graphite, iron Fe, magnetite FesC^). We 
then approximate the calculated absorption profiles Q a bs(A) with eq{T] For submicron-sized 
dust, as shown in Figured], while Q D « 1.6 is a reasonable estimation for both dielectric 
and metallic dust, the cut-off wavelength for a given size a is generally shorter for dielectric 
dust (A ~ naj^ than that for metallic dust (A Q ~ 2na). For MgS dust of size a, we 
thus adopt Q Q = 1.6 and A Q = ira. We note that, as can be seen in Figured] eqfj] with 



3 For amorphous silicate dust larger than ~ 0.07 /mi in size, A is much smaller than ira. This is due to 
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Q = 1.6 and A Q = ttcl (or A Q = 2na) almost always overestimates the actual UV/visible 
absorptivities of all six dust species, particularly at A < A Q . This suggests that the MgS 
temperatures derived using the UV/visible absorptivity of eqJT] would be overestimated and 
thus one would underestimate M Mg s (since it is most likely that eq{T] also overestimates the 
actual UV/visible absorptivity of MgS). 



In the IR, we will use the dielectric functions of Mg 9 Fe .iS measured by lBegemann et al. 



( 11994J ) in the 10-200 /im wavelength range to calculate the absorption efficiencies Q a b s (A) °f 



MgS. At longer wavelengths we will extrapolate assuming Q l ^ s (\) oc A~ 2 . We will consider 
both sphe rical dust and dust with a CDE (continuous distribution of ellipsoids) shape dis- 
tribution (IBohren &: Huffmanlll983l ). Finally, we smoothly join Q^ S (X) and Qabs(^) through 



Qabs(A) = (1 - £r) QZW + 6r<&(A), 912 A < A < 1 cm (2) 

^ir(A) = min [l, (A/lO/im) 3 ] . (3) 

The absorption efficiency synthesized from eqf2] is dominated by Q^sW at A < 5 jum and 
by Qj£,(A) at A > 10 fim (see Figure^. 



3. The Tester: HD 56126 

A successful candidate carrier should be able to explain the observed 30 /zm feature in all 
sources. A failure in a single source would be sufficient to rule out the candidate. To examine 
whether the carriers can account for the observed feature strength, we choose HD 56126, a 
prototypical source of the 21 fim and 30 /im features and one of the best studied post-AGB 
stars, as the tester. 



Van Winckel &: Reynierd (120001 ) performed a homogeneous photospheric abundance anal- 
ysis for HD 56126, and derived the abundances of sulfur and carbon to be S/H ph 4.07 x 10~ 6 , 
and C/H pa 4.47 x 10~ 4 . HD 56126 has a radius of r* pa 49.2 r Q (r Q is the solar radius), and 
a luminosity of pa 6054 L Q (L Q is the solar luminosity). We approximate the stellar ra- 
diation by the Kurucz model atmospheric spectrum with T e g = 7250 K and logg = 1.0. We 
adopt a distance of d ~ 2.4 kpc from Earth to HD 56126. 

The size and mass of the circumstellar envelope of HD 56126, which are important 
for determi n ing th e quantity of MgS dust available in the envelope, are still controversial. 



Hony et al.l (120031 ) suggested that the dust is confined to the envelope with radius between 



the sudden rise in absorption caused by the electronic transition at A ~ 0.1-0.2 /im (jKim fc Martin 



1995). 
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1.2"-2.6" based on their mid-IR imaging at 11.9 /im. Their detailed radiative transfer model 
derived a circumstellar envelo pe mass of Mh ~ 0.1 6-0.44 M^, depending on the assumed 
gas-to-dust ratio (~ 220-600). iMeixner et al.l (120041 ) found a more extended envelope (with 
radius between 1.2"- 7") based on their CO J = 1-0 line emission images. They derived a 
much smaller mass for the envelope (Mh ~ 0.059 M ). 



4. How Much MgS Dust Is Required? 

By assuming the envelope is optically thinQ we calculate the temperature T(r, a) of 
MgS dust of size a at a distance r from the central star from the energy balance between 
absorption and emission 

/r n2 /» r°° 
7T / ™ 2 Q abs (a,A)F A *dA= / ™ 2 Q abs (a, \)4nB x (T[r, a]) dX , (4) 
V/Ir/ Jo Jo 

where F£ is the stellar atmospheric flux which is approximated by a Kurucz model. In 
Figured we show the equilibrium temperatures of MgS dust of a = 0.1, 0.3 fim as a function 
of distance from the central starjfl We see that for a < 0.1 /iin T decreases from ~ 110 K at 
the inner most shell to ~ 50 K at the outer boundary, significantly lower than that adopted 
by iHony et al.l (120021 ). T = 150 K. Larger dust has a lower T (e.g. for dust of a = 0.3 /zm, T 
is lower than that of a = 0.1 /im by about ~ 10 K). 

The IR emission per unit mass of MgS dust of size a is 

^Mgs(A) = / [3Qabs(a, A)/4ap Mg s]4vr J B A (T[r,a]) dn(r) / dr Anr 2 dr , (5) 

where pMgS ~ 2.84 g cm -3 is the mass density of MgS dust, r in ps 1. 2" and r niltl ~ 7" ar e 
respectively the inner and outer radius of the dust shell of HD 56126 (IMeixner et al.l 120041 ). 
and dn(r)/dr is the MgS dust spatial density distribution. We will consider two spatial 
distribution functions: (1) dn/dr oc r~ 2 , which means the mass loss of the central star is 



4 If the envelope is not optically thin the conclusion of this work would be strengthened since when exposed 
to an attenuated starlight radiation, the MgS model would require more MgS dust to account for the same 
30 /im feature strength. 

5 For dust in the size range of 0.01 < a < 0.1 /im, T is insensitive to a since both Q^ s /a and Q l ^ s /a are 
independent of a. This is because for A > 0.3 /xm where there is significant stellar radiation, Q^bs oc A oc a, 
while in the IR the dust is in the Rayleigh regime and therefore Q^bs K a - 
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consta nt during the whole process^ and (2) dn/dr oc r 1 , which was shown by Hony et al. 



( 120031 ) to closely fit the observed SED and mid-IR images (with the dust confined to a narrow 



zone of 1.2"-2.6" from the star). 

The power output per unit mass in the 30 jum band is calculated by integrating FM g s(A) 
over the entire band (with the continuum subtracted) 

E 3 Z m = I F MgS (\)d\ . (6) 

J 30/xm band 

The MgS mass required to account for the observed 30 /xm emission in HD 56126 is 

A/f rec l — pobs / 77i30^m /»\ 
JW MgS ~~ • C/ 30|tm/ -^MgS i V ) 

where ffgff re2x l0 36 erg s 1 is the total power emitted from the 30 /iin feature of HD 56126 



(IHonv et all 120031 ). 

We calculate M^ s and tabulate the results in Tabled] We see that for dust with a 
size a > 0.1 /xm, M^ s increases with a because the dust temperature T decreases with a 
(see Figure^]). But for a < 0.1 /xm, M^ s is insensitive to a (since T is insensitive to a as 
long as a is not in the nano size range and the dust will not undergo stochastical heating by 
single photons). The required MgS mass is also affected by the adopted dust spatial density 
distribution: the dn/dr oc r _1 distribution requires more MgS (since it places more MgS 
grains in the cool, outer envelope region). But in general, is in the order of several 

tens of 1O~ 5 M exceeding what would be available in HD 56126 roughly by one order of 
magnitude (see §5]). 

The 30 /xm feature arising from spherical MgS dust is too sharp to be comparable with 
the 30 /xm featu re seen in C-rich objects. This is why a CDE shape distribution is often 



i nvoke d (e.g. see lHony et al.ll2002l ). Unfortunately, the simple formula of iBohren fc Huffman 
( 119831 ) for calculating the absorption cross sections of CDE dust is only valid in the Rayleigh 
regime (i.e. a « A). There is no precise way to derive the UV/visual absorption efficiencies of 
CDE dust (which are crucial to determine its heating rate and equilibrium temperature). We 
therefore simply adopt the equilibrium temperatures calculated for its spherical counterpart. 
The MgS mass required to account for the observed 30 /xm emission power is derived from 



fir 



The much more complicated function adopted bv lMeixner et alJ ([200J) to reproduce the observed dust 



IR emission spectral energy distribution (SED) and mid-IR images of HD 56126 is actually similar to dn/dr oc 



7 



This also a ppears unreasonably too high in view th at the total dust mass in the envelope is just 



7.4 x 1CT 4 M Q (jHonv et al.ll2003t iMeixner et alll2004h 



ust 
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eqsH M and [7] with Q^ s calculated from the CDE shape distribution (IBohren fc Huffman 
19831 ) . As shown in Tabled! the CDE model requires ~ 30% less MgS mass than the spherical 
model, but still needs much more than what could be available in HD 56126. 



5. How Much MgS Dust Would Be Available in HD 56126? 



Let be the MgS dust mass available in the circumstellar envelope of HD 56126. 

Let M env , Mh, and Md US t be the total gas, hydrogen and dust mass in the envelope. If we 
know the sulfur abundance [S/H]^, the fraction of sulfur depleted in MgS MgS , and the 
gas-to-dust ratio r g / d of HD 56126, we can estimate M^ s from 



M ava 



A*M g sM H [S/H], 



^MgS 



M H = M wy /lA , M env = r g/d M dust 



where /iM g s is the molecular weight of MgS, and the factor of 1.4 accounts for He whose 
abundance is ~ 10% of H. 



While both iHony et al.l (120031 ) and iMeixner et al.l (120041 ) derived the circumstellar dust 
mass to be M dust « 7.4 x 10~ 4 M o , they adopte d a very different ga s-to-dust ratio r g / d . 
Assuming that all C atoms are locked up in CO, IMeixner et al.l (120041 ) estimated M env ~ 
0.059 M & from the millimeter interferometry images of the CO J = 1-0 line. This corresponds 
to a gas-to-dust ratio of r g / d m 7 5, lower by a factor of ~ 3 than the typical gas-to-dust ratio 
°f r g/d ~ 200 for carbon stars (I Jural Il986l ). IHony et al.l (120031 ) argued for a much higher 
gas-to-dust ratio, with r g / d ~ 220-600. 

I f we assume that all S atoms ar e tied up in MgS (i.e. 0M g s = 1)' whh [S/H]^ « 4.07 x 
10~ 6 J Van Winckel fc Reynierlboooh and M dust « 7.4 x 10~ 4 M Q , the total amount of MgS 
mass available in HD 56126 would be M^ s » 9.60 x 10~ 6 , 2.60 x 10" 5 , and 7.16 x 10" 5 M & 
respectively for r g / d = 75, 220, and 600 (see Table[2]). Even with r g / d = 600, M^ s appreciably 
exceeds what is available in HD 56126! Therefore, MgS is unlikely the carrier of the 30 fim 
feature. 



Discussion 



While the gas-t o-dust ratio of r g / d = 600 adopted by IHony et al.l (120031 ) may be in the 
high end, r g / d ta 75 ( Meixner et al. 2004 ) is probably too low@ The problem probably lies 



8 With r g / d « 75 (corresponding to M cnv rj 0.059 M Q ), the star (whose mass is ~ 0.6 M ) would just have 
a ZAMS mass of M* w 0.66 M . Such a low-mass star, even with a metallicity [Fe/H] = —1 can not evolve 
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in the fact that iMeixner et al.l (120041 ) could have underestimated the total gas mass M env by 
assuming that all carbon atoms are in CO. As a ma t ter o f fact, a considerable fraction of 



the C atoms are in atomic carbon (CI): iKnapp et al.l (120001 ) observed the 609 /zm 3 Pi — > 3 P 



line of CI in the envelope of HD 56126 and estimated the CI to CO abundance ratio to 

ion of the C atoms should be C2, C3, and CN 



be [CI/COl « 


0.4. In 


(Bakker et al. 


1996: 1 



Let us assume that CO, CI, and amorphous carbon dust are the major sinks of C (i.e. 
we neglect C2, C3, CN and PAHs). The C budget in the envelope of HD 56126 is 



[C/HL = [C/H] co + [C/H] CI + [C/H] AC 



(9) 



wher e [C/H]^ « 4.47 x 10 4 is the stellar C abundance of H P 56126 (I Van Winckel fc Reyniers 
203), [C/H] co = M co / [mcoMh], [C/H] ci » 0.4 [C/H] co fcnapp et al.ll200oh . and [C/H]~= 
Mac/ [fJ'cMu], are respectively the C abundances ti ed up in CO, CI and amorph ous carbon, 



wher e M r ,n ~ 5.54 x 10~ 4 M m ( IMeixner et al.ll2004f ) and M AC « 7.4 x 1O" 4 M (IHonv et al. 



2003 



Meixner et al.ll2004j ) are respectively the masses of CO and amorphous carbon in the 
HD 56126 envelope, /ico an d He are the molecular (atomic) weights of CO and C, respec- 
tively. We estimate the total hydrogen mass to be 



1.4M C o/Aico + M A c//i 



c 



[c/hl 



0.20 M, 







(10) 



env ~ 
MgS 



0.28 M Q and a gas-to-dust ratio of r g / d ps 
= 1), we obtain M^ s « 4.56 x 1O" 5 M . 



This corresponds to a total envelope mass M< 
380. With all the sulfur atoms in MgS (i.e. 
This is much smaller than M^ s , the amount of MgS mass required to explain the power 
emitted at the 30 fim feature (see Table 1). 

We should note that the actual available MgS mass could be even lower since it is 
likely that some S atoms are tied up in gas molecules such as SiS and CS. In the C-rich 



( Glasseold 


1996; 


Millar et al. 


2001) 



As illustrated in Figured! Q Q = 1.6 and A = 7ra (adopted for MgS) also overestimate the 
UV/visible absorptivity of amorphous carbon. The mixture of amorphous carbon and MgS 
would have a lower temperature since the IR emissivity of amorphous carbon is much larger 
than that of MgS. Therefore one would expect a larger M^ s . We should also note that 



to the tip of RGB in the cosmic age (|Demarque et al.ll2008T ). not to mention the post-AGB phase. 

9 Goebel & Moseley (1985) argued that MgS could form via chemical surface reactions on carbon grains. 
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HD 56126 is not the strongest 30 /im feature source. The MgS budget problem would be 
more severe for objects with a stronger 30 /im band and c ooler MgS dust (e.g. IRAS 19454, 
a post-AGB star, requires T « 50 K; see lHony et al.ll2002l ). 



Zhukovska fc Gaill (120081 ) recently performed a detailed study of the condensation of 



MgS in the outflows from C-rich stars on the tip of the AGB. They found that MgS can 
only be formed by precipitating on pre-existing grains and therefore one would expect MgS 
to form as a mantle on the SiC core. However, these grains would exhibit a feature at 
~ 33-38 /jm which is not seen in C-rich objects. 



Hony et al.l (120031 ) and iHony Bouwmanl (120041 ) assumed an analytic formula for the 



UV/visible absorption efficiency of MgS dust, which differs from eq.l in that theirs have an 
appreciable absorbing power at 0.3 < A < 2 /mi: 




(11) 



With this kind of Q]^ s (\) and assuming a spatial distribution of dn/dr oc r" 1 (within a 
narrow zon e of 1. 2" -2. 6" from the star) for the dust and a gas-to-dust ratio of r g /d = 600, 
Hony et al.l (120031 ) performed detailed radiative transfer modeling calculations for HD 56126 
and found that the S abundance constraint was not violated. 



If we adopt the UV/visible absorption efficiency Q^ S (X) of IHony et al.l (120031 ). for MgS 
dust of a = 0.1 /im we would need M^ s « 8.38 x 10" 4 M© and 4.56 x 10~ 4 M respectively for 
the dust spatial distributions of dn/dr oc r _1 and dn/dr oc r~ 2 , far exceeding the maximum 
amount of MgS dust M^ s that could be available in HD 56126 (see Table 2)0 For a = 
0.01 /im, the required MgS mass is reduced to M^ s rj 8.01 x 1O~ 5 M and 4.69 x 10" 5 M Q 
respectively for dn/dr oc r~ l and dn/dr oc r -2 , and appears comparable to M^ s . However, 
we note that (1) M^| s is already an upper limit, and (2) we are not sure if eq.ll is a 
reasonable approximation for the UV/visible absorption properties of MgS dust of a = 
0.01 /im. As shown in Figure 1, the UV/visible absorption cut-off wavelength A D varies with 
dust size a. For the dust species listed in Figure 1, eq.l seems to be a more reasonable 
approximation than eq.ll for their Q^ S (X). We call on laboratory measurements of the 
UV/visible/near-IR dielectric functions of MgS dust. 



Finally, we calculate E^ s ( erg s 1 g 1 



the total power absorbed by one gram MgS 



10 Following Meixner et al. (2004), we assume the dust to be distributed in a broader zone, i.e., 1.2"-7" 
from the star. 
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dust of size a from 

Kt = ^ ln ^ out/r ; n) r o [3Q abs (a,A)/4ap Mg s] F^dX , ioidn/dr^r~ l ; (12) 

1 r out r in ^912 A 

r 2 1 Z" 00 

^abs = f / =[ 3 <3abs(M)/4ap Mg s] F x dX > for an/dr oc r~ 2 . (13) 

4 ri n r out J912 A 

If we assume the absorbed energy is a// radiated away through the 30 /iin feature, we would 
require a total MgS mass of = ^3o^ m /^abs — apparently, this is absolutely a lower 

limit since the absorbed energy will not be exclusively emitted from the 30 /im feature (a 
fraction of the energy will be radiated away at other wavelengths and through the continuum 
underneath the 30 /zm feature). In Table 3 we tabulate E^ s and M^g. We see in all cases 
Mjjgg > M^ s , indicating that we simply do not have enough MgS to account for the 
observed emission power. 

7. Summary 

We have investigated the hypothesis of MgS as a carrier of the prominent 30 /iin emission 
feature seen in numerous C-rich evolved objects, using HD 56126 as a test case. It is found 
that, in order to account for the enormous power emitted from this feature, one requires a 
much higher MgS dust mass than available in this object. We therefore argue that MgS is 
unlikely the carrier of the 30 /im feature. 

We thank S. Hony and A.K. Speck for very helpful discussions or comments. BWJ and 
KZ are supported by China 973 Program No. 2007CB815406. AL is supported in part by 
Spitzer Theory Programs and NSF grant AST 07-07866. 
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Fig. 1. — Absorption efficiencies Q a bs(A) of both dielectric (amorphous carbon, organic refractory, amor- 
phous silicate) dust and metallic dust (graphite, iron, magnetite) calculated from Mie theory for three sizes: 
a = 0.05 (red), 0.1 (green), 0.3 /im (blue). Also shown are the best fits (dashed lines) given by eq.l (with Q Q 
and A treated as free parameters). 
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Fig. 2. — Synthetic absorption efficiencies of MgS dust of a = 0.05 (red), O.f (blue), 0.3 /im (magenta). For 
a = 0.05/im we also show the UV/visible component Q^bs (g reen dot-dashed line) and the 1R component 
Oabs (black dashed line) alone. 

Table 1: M^ s (in unit of 1O _5 M ) — The MgS mass required to account for the power 
emitted in the 30 /im feature in HD 56126. 
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Fig. 3. — Equilibrium temperatures of MgS dust in the HD 56126 envelope. For dust with a < 0.1 /im, T 
is insensitive to a (see Footnote-[5]). 



Table 2: M^ s (in unit of 10 5 M Q ) — The MgS mass that could be available in the envelope 
of HD 56126, with the assumption that all S atoms are tied up in MgS. 
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Table 3: M^gs — The minimum amount of MgS dust required. This is obtained by assuming 
that the energy absorbed by a MgS grain is exclusively emitted through the 30 /im band (see 
§6). Also tabulated is E^ s (ergs _1 g~ 1 ), the power per unit mass absorbed by MgS dust. 
We see in all cases MSjS > MS,. 
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